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Abstract 

We study the interplay between coherent transport by tunnehng and diffusive 
transport through classically chaotic phase-space regions, as it is reflected in 
the Floquet spectrum of the periodically driven quartic double well. The tun- 
nel splittings in the semiclassical regime are determined with high numerical 
accuracy, and the association of the corresponding doublet states to either 
chaotic or regular regions of the classical phase space is quantified in terms of 
the overlap of the Husimi distribution with the chaotic layer along the sepa- 
ratrix. We find a strong correlation between both quantities. They show an 
increase by orders of magnitude as chaotic diffusion between the wells starts 
to dominate the classical dynamics. We discuss semiclassical explanations for 

this correlation. 
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I. INTRODUCTION 



In mesoscopic systems with a bistable potential, two fundamentally different modes of 
transport can occur simultaneously: Coherent quantal transport by tunneling, and deter- 
ministic diffusion through chaotic regions interweaving the classical phase space. There is a 
strong mutual influence of these two processes. While the classical phase-space structure is 
clearly reflected in the quantal spectrum [^]-|^, localization and delocalization due to quantal 
interference, in turn, effectively alter the boundary conditions for the quantal counterpart 
of classical phase-space diffusion [^|-|TD|. 

The idea that chaotic layers separating symmetry-related regular regions can give rise 
to quantal tunneling, just as potential barriers do, has been introduced by Davis and Heller 
0, who dubbed this notion "dynamical tunneling". The fingerprint of dynamical tunneling 
in the spectrum has been studied by Bohigas et al. P,^]. They utilized the occurrence of 
tunnel doublets with exponentially small splittings as a filter to separate "regular" quantal 
eigenstates (i.e., those localized in classically regular phase-space regions) from "chaotic" 
ones (in an analogous sense). While their work was based on a time-independent, two- 
dimensional nonlinear oscillator, the first inquiry into driven tunneling in a one-dimensional. 



bistable system was undertaken by Lin and Ballentine |^ and Grossmann et al. p!^-P 
and subsequently by Plata and Gomez Llorente |]T3 



Up to now, the only attempt towards an analytical understanding of tunneling between 
quantum states localized on symmetry-related tori (or chaotic regions) in multi-dimensional 
systems has been made by Wilkinson and Hannay ||T6| , p!7|| . Using semiclassical methods, 
they express the tunnel splitting as a function of the classical action along paths connecting 
the classical manifolds that support the doublet states. However, their theory does not 
provide a unified description spanning the transition from regular to chaotic states. It is 
this transition which particularly interests us here. 



Our present work (a preliminary account of which has been published in ref. 



IS 



similar in scope to that of Bohigas et al. mentioned above, but devoted to a different sys- 
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tern class. Applying torus or EBK (Einstein-Brillouin-Keller) quantization for periodically 
driven systems as well as the random-matrix theory for mixed (regular and chaotic) 
systems BJ20|] to the Floquet spectrum of the periodically driven double well, we expect the 
following scenario: The classical phase space of this system, for not too strong driving, is 
characterized by a pair of symmetry-related regular regions near the minima of the unper- 
turbed potential, embedded in a chaotic layer along the unperturbed separatrix. Quantal 
eigenstates localized in the regular regions form tunnel doublets with exponentially small 
splittings. These doublets will break up, and their splittings will reach values of the order 
of the mean level spacing, as soon as the associated pair of quantizing tori dissolves in the 
chaotic layer. Therefore, the tunnel splitting should depend strongly on the nature — regular 
or chaotic — of the classical phase-space region to which the corresponding pair of eigenstates 
predominantly belongs. In order to allow for a numerical test of this conclusion, we shall 
quantify this association in terms of the overlap of a suitable phase-space representation of 
the eigenstates with the chaotic layer. 

We emphasize that our hypothesis relies crucially on the specific phase-space topology 
of driven bistable systems. That is, the decisive property that prevents states within the 
chaotic layer from forming tunnel doublets is not a positive local Lyapunov exponent of the 
classical dynamics, but rather the possibility of phase-space transport across the symmetry 
plane. For example, in the converse case — symmetry-related chaotic islands embedded in 
regular regions — one expects tunnel-splitted level pairs which, in turn, repel each other in 
the same way as the single levels of the individual chaotic regions would do. Furthermore, 
it should be kept in mind that the disintegration of a classical torus is not an abrupt all-or- 
nothing event. Rather, this invariant manifold passes through a fractal, intermediate stage 
("cantorus" [BlI]) before it disappears, and even after that, a "shadow" remains in the form 



of a distinct, repelling structure within the chaotic sea ("vague torus" [Q). Due to the 
finite phase-space resolution of the quantal dynamics, this transition is smeared out even 
further in the parameter space of the corresponding quantum system. 

We give a brief review of our working model, its symmetries and classical dynamics, in 



Section |T[ In Section |TT|, we introduce the methods and key quantities used in our numerical 
study and interpret our results on basis of the simple picture just sketched. Section ^ 
contains the summary and puts our findings in a broader perspective. 



II. THE MODEL: SYMMETRIES AND CLASSICAL DYNAMICS 

The harmonically driven quartic double well is described by the Hamiltonian 

H{x,p;t) = Ho{x,p) +Hi{x;t), 

ffo(-.p) = ^-i-= + g^-', (1) 

Hi{x;t) = xS cos{ujt). 

With the dimensionless variables used, the only parameter controlling the unperturbed 
Hamiltonian Hq{x,p) is the barrier height D. It can also be interpreted as the (approx- 
imate) number of doublets with energies below the top of the barrier. Accordingly, the 
classical limit amounts to letting D ^ oo. The driving is characterized by its amplitude S 
and frequency uj. In all our numerical studies, we kept barrier height and driving frequency 
fixed at the values D = 8 and u = 0.95, respectively. 

Besides invariance under time-reversal, x ^ x, p ^ —p, t —t, the unperturbed 
system possesses the spatial refiection symmetry x — >■ —x, p ^ p, t ^ t. For a general 
periodic driving, this twofold symmetry is destroyed and replaced by the discrete time- 
translation invariance under t —>■ t + 2tt /uj. On the quantum-mechanical level, this enables 
the application of the Floquet formalism p3|-p^. For the specific time dependence of a 



harmonic driving, the symmetry f{t + n/uj) = — /(t) restores a similar situation as in the 



unperturbed case: The system is now invariant against the operation |12-14,27 



TT 

P : p — i> —p, X — i> —X, t t -\ , (2) 

which may be regarded as a generalized parity in the extended phase space spanned hj x, p, 
and phase, i.e., time t mod {2tt/uj). As in the unperturbed case, this enables to separate the 



eigenstates into an even and an odd subset. Moreover, as far as the dynamics allows, the 
generalized parity gives rise to tunnel doublets with the corresponding pair of eigenstates 



residing within the vortex tubes |19| formed by the motion within the potential wells of the 
unperturbed system. These tunnel doublets, in particular, are the subject of our study. 

The unperturbed system, as given by Hq{x,p), exhibits the generic phase space of a 
bistable system (Fig. la). Its structure is even independent, up to a rescaling of x and p, of 
the parameter D, provided D > 0. Apart from a constant of the order of unity, D gives the 
phase-space area enclosed by the two loops of the separatrix. It represents the characteristic 
action of the unperturbed system to be compared, in particular, with h. 

For the parameter values of the driving that are relevant for the present study, u = 
0.95 and < S < 0.2, the phase space of the unperturbed system is modified mainly 
by two additional features [^]: The onset of chaos in the vicinity of the separatrix, and 
the growth of the first nonlinear resonance zone (Figs. lb,c). The chaotic layer along the 
separatrix develops out of the homoclinic tangle, the intricate interweaving of the stable and 
the unstable manifolds originating at the hyperbolic fixed point near the top of the barrier 
PP| . It is present for all S > 0, but covers a substantial part of phase space only for S ^ 0.1. 

Nonlinear resonance arises when (an integer multiple k of) the period of the driving equals 
the period of the oscillation within each well at some energy. Each resonance generates a 
chain of k alternating elliptic and hyperbolic fixed points, embedded in the regular regions 
within each well. For example, the vortex tube surrounding, in the extended phase space, 
the elliptic fixed point of the first resonance, winds around the vortex tube corresponding to 
nonresonant motion once per period of the driving. Within the parameter interval 10^^ ^ 
S ^ 10^^, this first resonance zone grows from near invisibility to almost the full area of 
the regular region within each well. The higher resonances occupy only extremely narrow 
filaments in phase space which accumulate at the top of the barrier. Therefore, the border 
zone between the chaotic layer and the regular regions, formed by the higher resonances, is 



also very narrow, and in fact comes close to a sharply defined borderline As a further 
consequence, the chaotic layer "gains ground" in a quite smooth manner, as S increases. 



The intricate, self-similar layering of the border zone implied by the KAM theorem does 
not visibly affect this process. 

For larger values of the driving amplitude, higher-order features, such as secondary is- 
lands, appear and render the phase-space structure increasingly complicated (Fig. Id). As 
we shall see below, the growth of the first resonance leaves no significant trace in the tunnel 
splittings, while the consequences of the spreading of the chaotic layer are drastic. 

III. NUMERICAL METHODS AND RESULTS 

The basic ingredient required to study the quantum-mechanical aspects of a periodically 
driven system is the Floquet operator p3|- p^] , i.e., the unitary propagator that generates 
the time evolution over one period of the driving force, 

U = Texpl--J^ dtH{t)j, (3) 

where T denotes time ordering. Its eigenvectors and eigenphases, referred to as Floquet 
states and quasienergies, respectively, can be written in the form 

|^„(t))=e--*|0„(t)), (4) 

with 

\(l)a{t + 2n/iu)) = \Mt))- 
From a Fourier expansion of the | (pa {t) ) , 



-inuit 
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n 

Ca,n) = — / dt\(f)a{t))e 

In Jo 



(5) 



inuit 
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it is obvious that the quasienergies come in nu, n = 0,±1,±2,..., 

where each member corresponds to a physically equivalent solution. Therefore, all spectral 
information is contained in a single "Brillouin zone", —uj/2 < e < uj/2. 
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By inserting the eigenstates (4) into the Schrodinger equation, Fourier expanding, and 
using the representation in the eigenstates of the unperturbed Hamiltonian, i^ol^fc) = 



Ekl'^k), a matrix eigenvalue equation [p3|,p6 



is derived, where 



(6) 



n' k' 



Hn,k;n',k' = {Ek " nUj)Sn-n'Sk~k' 

r, Sn~l~n' + ^n+l-n' 
+ J Xk,k' ^ 5 

Cn,k = I C„), 

It is this matrix eigenvalue equation which we solve numerically to obtain the Floquet states 
and the quasienergies. 

The invariance of the system under the generalized parity discussed in the previous 
section (see Eq. (2)) is of considerable help in the treatment of Eq. (6). The eigenvalue 
equations for the subspaces spanned by the even and odd eigenvectors decouple completely, 
and the matrices to be diagonalized take the block structure 

/ ■ ... . \ 
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for the even (e) and the odd (o) subspace, respectively, where 
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2^0,1 xo,z 3;o,5 
2^2,1 2;2,3 2:2,5 
3^4,1 2:4,3 X4,5 



V 



Xo 



Xifi Xi,2 a:i,4 

3^3,0 2:3,2 3^3,4 
2^5,0 2:5 2 3:5 4 



V 



All the matrices involved are formally infinite and have to be truncated according to the 
required numerical accuracy. Evidently, exploiting the generalized parity reduces the linear 
dimensions of vectors and matrices by a factor of two. 

The tunnel doublets we are looking for are then composed of one state each from the 
even and the odd subspace, and the corresponding tunnel splittings are given by 



(8) 



Due to its periodic, Brillouin-zone structure, the quasienergy spectrum gives no hint how to 
order the doublets. Instead, it is possible to associate a mean energy to the eigenstates by 



the relation 26 



Eo. = ^ dt ( I H{t) I ), (9) 

/TT JO 

which is defined on the positive real axis and thus provides a basis for an ordering of the 
Floquet states. Using this order, we are able to pick out the doublets from the "ground- 
state" one upward to that immediately below the top of the barrier. In the following, we 
shall concentrate on this set of states, which for D = 8 are 18 in number. 

A quantity which facilitates making contact with the classical dynamics is the Husimi 
distribution pT|J5^ . It is defined as the overlap of a given state with a minimum-uncertainty 



(coherent) state localized at a position in phase space. 



with 



Q,(x,p;t) = -^|(e|^.(t))|', (10) 



oo 



n\ 



n=0 

where ^ = x + ip, and | n) is an eigenstate of the harmonic oscillator (p^ + x^)/2. The Husimi 
distribution represents a quantum state as a proper phase-space probability distribution with 
the highest resolution allowed by the uncertainty principle. 

In particular, the Husimi distribution can be used to quantify the notion of an eigenstate 
residing predominantly in a specific (regular or chaotic) region of classical phase space. We 
define the mean overlap of a state | ipait) ) with the chaotic layer as 

_ r'27v/uj poo poo 

Ta = — dt dx dpQa{x,p;t)T{x,p;t). (11) 

ZTT Jo J-oo J-oo 

Here, r{x,p;t) denotes the characteristic function for the chaotic region in the vicinity of 
the separatrix. Since the Husimi distribution forms a normalized probability distribution 
over phase space, we have < < 1. The characteristic function can be determined 
numerically, e.g., by letting a trajectory started anywhere in this chaotic region "tick" boxes 
in a coarse-grained phase space of the desired resolution. 

We start the discussion of our numerical results with a phenomenology of the Floquet 
states involved in driven tunneling. There are three aspects to be mentioned, the dependence 



on the quantum number (we choose the doublets 1 to 7), on the amphtude S of the driving, 
and on time t mod (271/00). For small enough S, e.g., S = 10~^, the Floquet states are 
nearly identical with the unperturbed eigenstates (Fig. 2). They are localized on tori with 
increasing characteristic action, embedded in the regular regions in the wells. The first 
qualitative changes to occur with S taking larger values are due to the growth of the first 
resonance (see the previous section). For S = 10^^ (Fig- 3), the doublet states with / = 2, 4, 5 
take shapes that indicate their localization on tori that belong to the first resonance rather 
than to nonresonant motion in the wells; the characteristics of torus quantization remain. 
At 5 = 0.2, the chaotic layer occupies an appreciable part of phase space. For the higher- 
lying doublets, the distorted toroidal shape of the Husimi distribution then gives way to a 
more rugged form which less closely resembles a smeared-out torus (Fig. 4). A glance at 
the time dependence within one period complements these observations (Fig. 5): For a state 
pertaining to the first resonance, e.g., the two maxima of the Husimi distribution rotate 
clockwise, with the frequency of the driving, around the respective potential minima, as do 
the corresponding classical vortex tubes. 

The central result of our study is presented in Fig. 6, where we compare the S dependence 
of the tunnel splittings (part a) with that of the overlaps F/ e (part b; we could quite as 
well have chosen F^^o)- Looking at the tunnel splittings first, a relatively clear-cut picture 
emerges: For 5* ^ 10^^, the tunnel splittings do not deviate significantly from their unper- 
turbed values (apart from an irregularity in A2 which roughly coincides with a region of rapid 
growth of the first resonance). They increase roughly exponentially with the quantum num- 
ber /, from a value Ai = 5.20 x 10~^^, in good agreement with the corresponding semiclassical 



(instanton-method, cf. ref. H) estimate H Af = Jl28D/n exp(-16D/3) = 5.33x10 



-18 



to Ag = 5.59 X 10~^. For S ^ 10^'^, the splittings start one by one to grow exponentially 
with S, from the ground-state doublet (/ = 1) upwards, so that at S* = 0.2, the range of 
the splittings has shrunk from 14 to three orders of magnitude. The 5* dependence of the 
overlaps shows the same qualitative features as that of the tunnel splittings: There is no sig- 
nificant deviation from the unperturbed values for S ^ 10~^, while for larger S, exponential 
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growth sets in from bottom to top of the doublet ladder, so that the F/^e successively join 
in an approximately single line. 

The qualitative agreement between the two respective groups of functional dependences 
indicates that there exists a strong correlation between the tunnel splittings and the overlaps 
with the chaotic layer. Furthermore, the steep exponential increase occurring in both quan- 
tities coincides with the onset of chaotic motion in the classical dynamics, whereas, e.g., the 
transfer of phase-space area from nonresonant motion to the first resonance hardly leaves 
any trace in the tunnel splittings. Insofar, the simple picture sketched in the Introduction 
is confirmed. Details of our expectation, however, need to be revised. 

In particular, the notion that each splitting widens up individually as the corresponding 
quantizing torus resolves, is not unambiguously corroborated by the data. It would imply 
that the transitions to a large splitting occur from top to bottom, i.e., first for the doublet 
localized on the outermost torus, the one with the highest mean energy. Indeed, if this 
transition is assessed from the splittings passing a certain absolute threshold, say A; = 10~^, 
that order is roughly followed (at least for those doublets we can keep track of that far). If, 
however, the point of onset of exponential growth, visible in a logarithmic plot, is taken as 
the criterion, the order is reversed. 

Another remarkable fact is that the widening of the splittings, and the concomitant 
change in character of the eigenf unctions, as a function of S, is a continuous process that 
can only vaguely be associated with the decay of a KAM torus, taken as as a discrete event. 
Even doublet states overlapping by 70% with the chaotic layer may still show a relatively 
small splitting and exhibit the signature of a regular state in their spatial structure and time 
dependence (see, e.g.. Figs. 4 and 7). It remains to be clarified whether this retarded decay 
of the tunnel doublets corresponds to the gradual disintegration of classical tori via cantori 
and vague tori. 

Unfortunately, our numerical means do not allow to extend the data beyond S — 0.2. 
Therefore, the eventual saturation of the growth of the tunnel splittings, when they reach 
the order of magnitude of the mean level spacing, could not be studied. Furthermore, even 
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if the double well with a scaled barrier height D — 8 is the closest approach to the classical 
limit we could practically afford, the quantal dynamics of this system is not yet sensitive 
enough for the details of the classical phase space to allow for a conclusive numerical test of 
a semiclassical description. 

IV. SUMMARY 

In which way does the onset of classical chaos in the vicinity of the separatrix influence 
the tunnel splittings in a periodically driven bistable system? The random-matrix theory 
for mixed (chaotic and regular) systems, together with semiclassical considerations, suggest 
the following simple answer: A tunnel splitting undergoes a transition from exponentially 
small values to a size comparable with the mean level spacing, as the corresponding pair 
of symmetry-related, quantizing tori resolve in a common chaotic region. In order to test 
this hypothesis, we quantifled the association of an eigenstate to the chaotic layer along 
the separatrix in terms of the overlap of its Husimi representation with this phase-space 
region. We found a striking qualitative agreement between the functional dependences of 
the splittings and the corresponding overlaps on the amplitude of the driving force. Both 
groups of quantities show transitions to steep exponential growth as the chaotic layer along 
the separatrix spreads at the expense of the regular regions around the potential wells. The 
above hypothesis does, however, imply more, namely that the transitions to a large splitting 
for the individual doublets occur in a specific order given by the fate of the corresponding 
tori. For the driven double well this order is roughly from top to bottom on the energy axis 
(where "energy" should be read as "mean energy"). This is indeed what we observe if the 
widening of the splittings is assessed from their growing beyond some absolute threshold 
value. If, however, the onset of exponential growth is taken as the criterion, the order turns 
out to be from bottom to top. As a further unexpected detail, the regular character of 
doublet states survives until far into a parameter region where they are already located 
amidst the chaotic layer. 
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On first sight it may appear strange that a static quantity hke the overlap with the chaotic 
layer should bear on a transport phenomenon such as tunneling. Prom a semiclassical point 
of view, this may become understandable if one assumes that a path-integral expression for 
the tunnel splitting becomes dominated by contributions from paths which pass through the 
chaotic layer, as soon as chaotic diffusion enables significant classical transport between the 
wells. Our study, which is intended solely as a first numerical survey, should encourage to 
cast this notion in a precise semiclassical form. Related topics which deserve further investi- 
gation are the time-domain aspects of tunneling through the chaotic layer, the occurrence of 
crossings between tunnel doublets and their influence on the dynamics. Another interesting 
open question is whether states which reside within the chaotic layer but appear regular in 
character, can be associated with vague tori of the classical dynamics. 



13 



REFERENCES 

[1] G. M. Zaslavsky, Sov. Phys.-Usp. 22, 788 (1979) [G. M. Zaslavsky, Usp. Fiz. Nauk 129, 
211 (1979)]. 

[2] O. Bohigas and M.-J. Giannoni, Lecture Notes in Physics vol.209 (Springer, Berlin, 
1984), p. 1. 

[3] M. V. Berry, Proc. R. Soc. A 400, 229 (1985). 

[4] G. Casati, B. V. Chirikov, F. M. Izrailev, and J. Ford, Lecture Notes in Physics vol.93 
(Springer, Berlin, 1979), p. 334. 

[5] S. Fishman, R. E. Prange, and D. R. Grempel, Phys. Rev. Lett. 42, 509 (1982); Phys. 
Rev. A 29, 1639 (1984). 

[6] D. L. Shepelyansky, Phys. Rev. Lett. 56, 677 (1986). 

[7] M. J. Davis and E. J. Heller, J. Chem. Phys. 75, 246 (1986). 

[8] O. Bohigas, S. Tomsovic, and D. UUmo, Phys. Rev. Lett. 64, 1479 (1990); Phys. Rev. 
Lett. 65, 5 (1990). 

[9] S. Tomsovic and D. UUmo, Tunneling in the Presence of Chaos, preprint, 1990. 

[10] O. Bohigas, S. Tomsovic, and D. UUmo, Phys. Rep. 223, 43 (1993). 

[11] W. A. Lin and L. E. Ballentine, Phys. Rev. Lett. 65, 2927 (1990); Phys. Rev. A 45, 
3637 (1992). 

[12] F. Grossmann, P. Jung, T. Dittrich, and P. Hanggi, Z. Phys. B 84, 315 (1991); Phys. 
Rev. Lett. 67, 516 (1991). 

[13] F. Grossmann, T. Dittrich, P. Jung, and P. Hanggi, J. Stat. Phys. 70, 229 (1993). 

[14] F. Grossmann, T. Dittrich, and P. Hanggi, Physica B 175, 293 (1991). 

[15] J. Plata and J. M. Gomez Llorente, J. Phys. A 25, L303 (1992). 

14 



[16] M. Wilkinson, Physica 21D, 341 (1986); J. Phys. A 20, 635 (1987). 
[17] M. Wilkinson and J. H. Hannay, Physica 27D, 201 (1987). 

[18] R. Utermann, T. Dittrich, and P. Hanggi, in Proceedings of the XXth International Con- 
ference on Low Temperature Physics, Eugene, 1993, edited by R. J. Donnelly (North- 
Holland - Elsevier, Amsterdam, 1993). 

[19] H. P. Breuer and M. Holthaus, Ann. Phys. (New York) 211, 249 (1991). 

[20] M. V. Berry and M. Robnik, J. Phys. A 17, 2413 (1984). 

[21] L. E. Reichl, in The Transition to Chaos: In Conservative and Classical Systems: Quan- 
tum Manifestations (Springer, New York, 1992), Chaps. 3.9 and 9.5.1, and refs. therein. 

[22] W. P. Reinhardt, J. Phys. Chem. 86, 2158 (1982); R. B. Shirts and W. P. Reinhardt, 
J. Chem. Phys. 77, 5204 (1982). 

[23] J. H. Shirley Phys. Rev. 138B, 979 (1965). 

[24] H. Sambe, Phys. Rev. A 7, 2203 (1973). 

[25] N. L. Manakov, V. D. Ovsiannikov, and L. P. Rapoport, Phys. Rep. 141, 319 (1986). 
[26] S. Chu, Adv. Chem. Phys. 73, 739 (1989). 
[27] A. Peres, Phys. Rev. Lett. 67, 158 (1991). 

[28] L. E. Reichl and W. M. Zheng, in Directions in Chaos, vol.1, edited by H. B. Lin (World 
Scientific, Singapore, 1987), p. 17. 

[29] A. J. Lichtenberg and M. A. Lieberman, Regular and Stochastic Motion, Vol. 38 of Appl. 
Math. Sci. (Springer, New York, 1983). 

[30] H. P. Breuer, K. Dietz, and M. Holthaus, Physica D 46, 317 (1990). 

[31] K. Husimi, Proc. Phys. Math. Soc. Jap. 22, 264 (1940). 



15 



[32] R. J. Glauber, in Quantum Optics, edited by A. Maitland (Academic, London, 1970). 

[33] S. Coleman, in The Whys of Subnudear Physics, Proceedings of the International School 
of Subnuclear Physics, Erice, 1977, edited by A. Zichichi (Plenum, New York, 1979). 



16 



FIGURES 

FIG. 1. Classical phase-space portraits of the periodically driven double well at phase of the 

driving, for various values of the driving amplitude. The parameter values are -D = 8, u; = 0.95, 
and (a) S = 0,{h)S = 10'^, (c) S = 0.02, (d) S = 0.2. 

FIG. 2. Contour plots of the Husimi distributions for the Floquet states |'02,e(O)) (a) and 
|V7,e(0)) (b), at 5 = 10-5. 

FIG. 3. Contour plot of the Husimi distribution for the Floquet state | ■05,e(O) ) at S = 0.02. 
It should be compared with the corresponding classical phase-space portrait. Fig. Ic. 

FIG. 4. Contour plot of the Husimi distribution for the Floquet state | '07,e(O) ) at S = 0.2. 

FIG. 5. Contour plot of the Husimi distribution for the Floquet state | ■05,e(i) ) at S = 0.02 
and (a) cot = 7r/4 and (b) cot = it/ 2 (b), compared to the corresponding classical phase-space 
portraits at the same parameter values and phases (parts (c) and (d) , respectively) . For the initial 
states {ut = 0), see Figs. Ic and 3, respectively. The states at later times ut = mr/4 are related to 
those shown by simple phase-space symmetries, as implied by the generalized parity, Eq. (2). 

FIG. 6. Tunnel splittings (a) and overlaps with the chaotic layer (b) for the seven lowest 
tunnel doublets, as functions of the amplitude of the driving. 

FIG. 7. Contour plot of the Husimi distribution for the Floquet state | ip7,e{t) ) at S = 0.02 
and (a) = 7r/4 and (b) ut = 7r/2. For the initial state, see Fig. 4. The states at later times 
cot = mr/A are related by simple phase-space symmetries to those shown here, as implied by the 
generalized parity,Eq. (2). 
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